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Abstract 

The analytical solutions reported in our previous paper are given as series of hyperge- 
ometric or Coulomb wave functions. By using them, we can get the Teukolsky func- 
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tions analytically in a desired accuracy. For the computation, the deep understanding of 
their properties is necessary. We summarize the main result: The relative normalization 
between the solutions with a spin weight s and — s is given analytically by using the 
Teukolsky-Starobinsky (T-S) identities. By examining the asymptotic behaviors of our 
solution and combined with the T-S identities and the Wronskian, we found nontrivial 
identities between the sums of coefficients of the series. These identities will serve to 
make various expression in simpler forms and also become a powerful tool to test the 
accuracy of the computation. As an application, we investigated the absorption rate and 
the evaporation rate of black hole and obtain interesting analytic results. 
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1 Introduction 



In our previous paper [1], we reported the analytic solutions of Teukolsky equation [2] 
which consist of two type of series; one is given in the form of series of hypergeometric 
functions (hereafter we call it the hypergeometric type solution) and the other is given 
in the form of series of Coulomb wave functions (hereafter we call it the Coulomb type 
solution which was first given by Leaver [3]). The hypergeometric type solutions are shown 
to be convergent in the region except infinity with all finite e = 2Mu, M being the mass 
of black hole and uj an angular frequency. The Coulomb type solutions are convergent in 
the region far above the outer horizon for all finite e. We showed that the matching of 
these two types of solutions is perfect in the intermediate region where both solutions are 
convergent. In this way, we found the solutions which are convergent in the entire region 
of r for all finite e. In addition, we showed that our solutions are suitable to get the 
Teukolsky functions in the e expansion or the numerical computation. We also presented 
analytical solutions of Regge- Wheeler equation[4]. 

We showedfl] that the solution can be written in various forms so that we have to 
choose some special ones. Furthermore, it is necessary to investigate further the properties 
of solutions to calculate the Teukolsky functions, because the solutions are given as series. 
For this purpose, we examined the asymptotic behaviors of the solutions and compared 
with the general result derived before. We made a systematic study on the incoming 
solution on the outer horizon. 

As for the incoming solution on the horizon, we take the same hypergeometric type 
one as in Ref.l, for which the boundary condition on the horizon is trivial. For the 
outgoing solution, we take a different form. As for the Coulomb type solutions, we take 
different forms from that given in Ref.l. In this paper, we choose a variable z = uj{r — r + ) 
which is different from the definition in Ref.l where we tookz = uj(r — r_). Accordingly, 
the solutions adopted here are different. We took them because z reduces to oor in the 
Schwarzschild limit and thus the solutions for large r can be compared with the previous 
analysis by Tagoshi and Nakamura, and Sasaki and companies [5] for Schwarzshild case 
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and to our solutions of the Regge- Wheeler equation[4]. Therefore, some of solutions given 
here are new ones. 

In order to investigate the properties of the solutions, we mainly examine the Teukolsky- 
Starobinsky (T-S) identities [6], [7] and the Wronskian[6]. Although the solutions of Teukol- 
sky equation satisfies the T-S identities automatically, it is non trivial to see directly that 
our solutions satisfy them because our solutions are given as series. Fortunately, we were 
able to show analytically that the incoming solution satisfies one of T-S identities. By 
using this, we can fix analytically the relative normalization between solutions with a spin 
weight s and — s. It was hard to prove that our solution satisfies another T-S identities 
directly. In order to investigate the T-S identities and the Wronskian, we considered the 
asymptotic behaviors. For this, we derived the expressions of asymptotic amplitudes on 
the horizon and at infinity from our solution. By examining the T-S identities and the 
Wronskian, we found two very nontrivial identities among sums of coefficients. Since 
the solutions are defined as series, various quantities expressed by Teukolsky functions 
necessary involve the sums of coefficients which are calculated by solving the three term 
recurrence relation. These identities will serve to make the expressions of them in com- 
pact forms form which we may be able to have a direct incite of the physical meaning. 
As an example, we considered the absorption rate and the evaporation rate of black hole 
and obtained compact forms of them. For computing Teukolsky functions, we will make 
either the e expansion or the numerical computation of coefficients. In order to test the 
accuracy of the computation, we can use these identities. Since the sum of coefficients in 
these identities are proportional to the asymptotic amplitudes on horizon or at infinity, 
the check of these identities gives a direct check of accuracy of Teukolsky functions. 

With these understandings of the solutions given in this paper, it is now in the position 
to make the numerical study of Teukolsky functions. By solving the three term recurrence 
relation, we can in principle obtain Teukolsky functions with a desired accuracy. The 
further numerical study is in pressing need. 

In Sec. 2, we give the outline of solutions of Teukolsky equation and new forms of 
them. In Sec. 3, various properties of coefficients the proportionality factor between the 
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hypergeometric type and Coulomb type solutions are given. The asymptotic behaviors 
on the outer horizon and at infinity are expressed directly by using our solution in Sec. 4. 
In Sec. 5, we discuss properties of our solution by using the T-S identities. By using the 
Wronskian, we discuss the conserved quantity and the energy conservation in Sec. 6. By 
using the T-S identities and the Wronskian, we derive two identities among the sums 
of coefficients in Sec. 7. The absorption rate and the evaporation rate by black hole are 
discussed in Sec. 8. Discussions are given in Sec. 9. The outlines of proofs of various 
formulas in the text are given in Appendix A. 



2 Analytic solutions 

The hypergeometric type solution which we discuss here is the same as the one given in 
Ref.l. However, we take different solutions for the Coulomb type from the one in Ref.l. 
We do not present the derivations and the proofs of various equations and statements 
which one can find in Ref.l or one can prove following the discussions given there. 

(a) Notation 

We start from the Teukolsky function [2], [6] 

T s ~ e-^e im tSr(0)sKim(r), (1) 

where s S] n (9) is a spin- weighted spheroidal function with a spin- weight parameter s and 
I is an angular momentum which satisfies / > max(| m |, | s |), and s Ruim( r ) is a radial 
function which we write R s for simplicity. Solutions of the radial equation are expressed 
as functions of variables x and z which are defined by 

x = (r — r+), z = u>(r — r_) = e/c(l — x) — enx, (2) 

€K 

where 

e = 2Mu, K=y/l- q\ q = ^ . (3) 

Here, M and a are the mass and the angular momentum of the Kerr black hole, re- 
spectively. The parameters r± are poles of A = r 2 + a 2 — 2Mr and is expressed by 

3 



r± = M ± \/m 2 — a 2 = e(l ± k)/2uj where r± give positions of the outer and the inner 

horizon. For later use, it is convenient to define 

e — mq , ma n k ( cu\ , 

t = -, k = cu = 2 - e+ 4 

k ' 2Mr+ 1 + k V e / + K J 



and 



e ± ^. (5) 



(b) Coefficients 

Solutions of the radial part are expressed in the form of series of hypergeometric 
functions and also in the form of series of Coulomb wave functions. Coefficients of these 
series satisfy the same three term recurrence relation [1], [8] 



<(s)< +1 (s) + (%{ 8 )a»{ 8 ) + tM<Ci(*) = 0, (6) 



where 



v( . ien(n + v + 1 + s + ie){n + u + 1 + s — ie){n + v + 1 + it) 

{n + v+ l)(2n + 2z/ + 3) ' 
f3 v n {s) = -A s -s(s + l) + (n + z/)(n + z/ + l) + e 2 + e(e-mg) 
e(e — mq)(s 2 + e 2 ) 
(n + v){n + v + 1) ' 
ie/t(n + z/ — s + ie)(n + v — s — ie)(n + v — it) 
ln{3} = (n + u)(2n + 2u-l) ' ( ] 

where \ s = E(s) — 2amu + a 2 uj 2 — s(s + 1) with E(s) being the eigenvalue of the spin- 
weighted function s Sj n (9). 

The series of coefficients converges if the transcendental equation for v 

R n {v)L n _ x {v) = 1, (8) 

is satisfied, where R n (u) and L n {v) are the continued fractions defined by 



T (jA n{ ~ > — (q) 
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The solutions are characterized v which we called the renormalized (shifted) angular 
momentum because v — I + 0(e 2 ) as shown in Ref.l. We can prove that if v satisfies 
Eq.(8), then —v — 1 also satisfies, i.e., R n (— v — l)£n-i( — v — 1) = 1. In other words, if 
v is a solution, then —v — 1 is also a solution. Note that E(s) is an even function of s, 
E(—s) = E(s) as shown by Press and Teukolsky[9]. By using this, we can prove that v 
is an even function of s 

v{s) = v(-s), (10) 

which is important and enables us to relate the solution of spin weight s to that — s by the 
Teukolsky-Starobinsky identities. This can be proved by noticing that the transcendental 
equation (8) contains f3^(s) and a combination ct^(s)^ +1 (s) which are even functions of 
s. Thus, the equation (8) is invariant under the change of s to — s and its eigenvalue v(s) 
is an even function of s. 

For v which satisfies Eq.(10), we find 

,• a n( S ) ,• a n( S ) ieK /-.is 

which enables to establish the region of convergence of solutions. 

(c) Analytic solutions 

(c-1) Hypergeometric type solutions 
The incoming solution is given by 

oo 

iC; S = A^i-xY^il-xt- £ <( S ) 

n=— oo 

xF(n + v + 1 — it, —n — v — ir; 1 — s — 2ie + ; x) , (12) 

where a£(s)'s are coefficients which satisfy the three term recurrence relation in Eq.(6) 
and A s is the normalization constant which we take for s > with the choice A- s = 1 
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A-s — 1, A s — C* [ 



/lu\ 2s r(l + s-2ie H 



r(z/+ 1 -s + ie) 



r i/+ 1 + s + ie) 



(13) 



\e/t/ T(l — s — 2ie„ 

where C s 's are Starobinsky constants (defined in Eqs.(49)-(52)). The relative normaliza- 
tion is determined so that the Teukolsky-Starobinsky identities are satisfied as we shall 



see later. It may be worthwhile to mention that R^~ l = R\ n . s so that the change v into 
—v — \ does not lead to a new solution. 

Another independent solution is the outgoing solution on the outer horizon and is 
given by 

Kut; S = A- S (iC;-J* (14) 

which is different from the one given in Ref.l. Solutions, R" n . s and i?g Ut . s form a pair of 
independent solutions of Teukolsky equation. 

By using large |n| behaviors of coefficients in Eq.(ll) and hypergeometric functions, 
we can prove (see Ref.l) that the solution in Eq.(12) is convergent for the entire complex 
plain except infinity for all finite e. If we restrict x to the physical region, the convergence 
region is < (— x) < oo, i.e., r < oo. 

The solution can be analytically continued by 



T3U A IEK 



K, s + Ro,s , (is) 



where 



R u . s = e- leK %x) u+it +{x-iy s - it+ 

r(l-s-2ie+)r(2n + 2z/+l) 



oo 



n ±^ r(n + v + 1 - ir)T{n + v + 1 - s - it) n[ } 

xx n F(—n — v — it, —n — v — s — ie; —2n — 2u; — ) . 

x 



(16) 



Then, Rq. s and R .^ form a pair of independent solutions. 

(c-2) Coulomb type solutions 
The solution is given by 

BT - z-i-d-tHY"** V ( iY {u+1 + S ~ te)n ^(c)F (z) 

R C ;s ~ Z (1 J J^J t) {l/+1 _ s + ie)n a n{s)F n+U!S {z), 



(17) 
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where (a) n = T(n + a)/T(a) and 

p ( z ) — e- iz (2z) n+ »z rin + iy + 1 ~ S + ie) 
t n+vA z) - e [Zz) z r(2n + 2z/ + 2) 

x$(n + v + l- s + ie, 2n + 2u + 2;2iz) . (18) 

This solution and the definition of z are different from those given in Ref.l. It is amusing 
to observe that coefficients appearing in the hypergeometric type solution are the same 
as those in the Coulomb type solutions[l],[9]. Thus, the renormalized (shifted) angular 
momentum takes the same value both for the hypergeometric type and the Coulomb 
type solutions. This fact is quite important for matching of these two different types of 
solutions in the region where both solutions are convergent. 

Another independent solution is obtained by changing v to —v — 1, 

Rd7" ■ (is) 

Solutions Rq. s and -R^ -1 form a pair of independent solutions of Teukolsky equation. 

Coulomb type solutions given above contain both the incoming and the outgoing 
solutions at infinity. Thus, solutions are decomposed into another pair of solutions, the 
incoming solution at infinity R v + . s and the outgoing solution at infinity R v _. s . Explicitly, 
we have 

R^ S = R^ S + R%, (20) 



where 



R» = re -^jAv+l-s) T ^+ l S _±^> e -iz z »+ie+, z _ eK ys-ie + 

oo 

x S i n a v n {s){2z) n ^{n + v+ 1 - s + ie,2n + 2v + 2;2iz), 



(21) 



R u _. s = 2 u e- we e" l ^ u+1+s) e lz z u+le +(z - €K)- s - it+ ]T i n 

n=—oo 

x ^ + 1 + S ~ ie \ n a u n (s)(2z) n V(n + v + 1 + s - ie, 2n + 2v + 2; -2iz) . 
\y + 1 — s + ie) n 



(22) 
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By using the large | n | behaviors of coefficients in Eq.(f f) and Coulomb wave func- 
tions, we can prove (see Ref.l) that solutions are convergent in the complex region of 
| z \> en for all finite e. If we confine to the physical region, the convergence region is 
r > r+. 

(c-3) Matching of two types of solutions 

Let us compare Rq. s and Rq. s . We first observe that they are solutions of Teukolsky 
equation and behaves like x v multiplied by a single- valued function for large x (z — enx). 
Thus, Rq. s should be proportional to R u c . s . The same is true for Rq.^' 1 and -R^" 1 . It 
should be noted that v — {—v — 1) = 2v + 1 is not an integer so that Rq. s and -R^ -1 are 
solutions with different characteristic exponents. 

We define 

Rl s = K v {s)R^ s . (23) 

Then, by comparing each power of x in the region 1 « x < oo where both solutions 
converge, we find 

{2eny u+s ~ f 2~ s i f T{l - s - 2ie + )V(f + 2u + l)Y(f + 2u + 2) 



K v {s) 



x- 



T(v + 1 - it)Y(v + 1 - s - ie)Y(f + v + 1 - s + ie) 
T(u + 1 + ir)r(iy + 1 + s + ie) 
r(f + v + 1 + ir)r(f + v + 1 + s + ie) 

/g, + i + .- fc ) y 

1 ^ (n-f)! 1 - s + ie) nK ' ' 



\n=r 



x f f ^ ( y+1 + 3 -' e )" ^( J )V 1 

' ^ (f-n)!(f + 2u + 2) n (u+l-s + ie) n nK ') 



\n=—oo 



(24) 

where f can be any integer and K„(s) is independent of the choice of f. 

(c-4) Solutions valid in the entire plane of r 

Incoming solution on the horizon 

We take the hypergeometric type expression for R? given in Eq.(12) for the incom- 
ing solution on the outer horizon which is convergent except infinity. By the matching 



between the hypergeometric type solution and the Coulomb type solution in Eq.(23), we 
have the Coulomb type expression for R\ n . s 

R^ s = A s e™[K„(s)R» c . )S + K^^R^ 1 } , (25) 

which is convergent in the region r > r + . The solution which is defined in this way is 
convergent in the entire region of r for all finite e. That is, we use the hypergeometric 
type expression in Eq.(12) around the outer horizon and the Coulomb type expression in 
Eq.(25) around infinity. In the intermediate region, we use either the expression by using 
Rq. s in Eq.(15) or the Coulomb type expression depending on the situations. 

Outgoing solution on the horizon 

By using R\ n ._ s1 the outgoing solution on the outer horizon is given by R^ nt ._ s = 
A~ s (Ry n ._ s )* . Thus, this solution is valid in the entire region of r for all finite e. The 
relative normalization between solutions with spin weight s and — s is fixed automatically 
following the normalization of R\ n . s - 

Upgoing solution 

The upgoing solution is the one which satisfies the outgoing boundary condition at 
infinity so that we can write 

K P;s = B S R-; S , (26) 

where S s 's are normalization constants and the relative normalization between B s and 
B_ s is fixed for s > with the choice B_ s = 1 

B. a = l, B s = C> 2s , (27) 

where C s 's are Starobinsky constants (defined in Eq.(49)-(53)). By using the relations 

R - v -\ = „. c -j^ sin7r(i/- s + ie) RV 

+ ' s sin7r(z/ + s - ie) +;s ' 

iqr 1 = ie^id.,, (28) 



we find 



t>v f> ' I . smir(v + s-ie) x e ™ sin ir(v - s + ie) , 



Then, if needed, we can express R^. s in terms of R\ n . s and R^t-s an d then we obtain the 
hypergeometric type expression. Now, we obtained the upgoing solution which is written 
by the hypergeometric type expression which is convergent in the region r < oo and the 
Coulomb type expression which is convergent for r > r + . Thus, we find the upgoing 
solution which is convergent in the entire region of r for all finite e. 



3 Some properties of coefficients and K u {s) 

In this section, we present various properties of coefficients and some relations between 
K u (s)'s which become important in the following discussion. Hereafter, we consider s = 
(masless scalar), s = ±1/2 (massless fermion, neutrino), s = ±1 (photon) and s = ±3/2 
(massless spin 3/2 particle) and s = ±2 (graviton). 
(a) Properties of coefficients 

In below, we take initial values of coefficients for s = 0, ±1/2, ±1, ±3/2, ±2 



(30) 



Then, from the three term recurrence relation, we have 



a u n {s) . 



(31) 



In addition, by using the explicit forms of a£(s), Pn( s ) an< ^ 7«( s ) an d the three term 
recurrence relation, we obtain the following relations: 

2 



<(-s) 



[v+l + s + ie) 



[v + l-s + ie) 



<0), 



(32) 



(33) 



<(-*Y 



■ly 



(u +1 + ir), 
{y + 1 - it), 



{y + 1 + s + ie) n 
{y + 1 - s + ie) n 



<{*) 



(34) 
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To prove Eq.(32), we see that <(-s) = \(n + u + l - s + ie)/(n + v + 1 + s + ie)\ 2 a v n (s) , 
Pn(~ s ) = Pn( s ) and ln(~ s ) = \(n + u + s + it) / (n + v - s + ie)\ 2 ^(s) and then, it is 
clear that a u n (—s) defined in Eq.(32) satisfies the three term recurrence relation for — s. 
Similar considerations will lead to Eqs.(33) and (34). 

(b) Useful relations for K u (s) 
The first one is 

K v (-s) 



K- v -i{-s) 

i r(i + s 



2ie, 



(en) 2s r(l - s - 2ie+) 



T(u+ 1-s + ie) 



T(v+ 1 + s + ie) 



2s 



(35) 



which can be proved directly by using the relation between a^(s) and a v n {— s) in Eq.(32). 
The second one is rather nontrivial, but can be proved directly for f = 0, 

2e«r(l + s + 2ie+)r(l - s - 2ie + ) sin7r(z/ - it) 



K v {s)K. u ^(-sY 



7T 



X 



sin7r(z/ — s + ie) 



sin 1-kv 



(36) 



where we used Eqs.(31) and (32). Since K„(s) is independent on the choice of f, the 
relation in Eq.(36) is valid for all integer value of f. These are useful relations which can 
be used to prove that the solution satisfies the T-S identity directly and to discuss the 
absorption rate and the evaporation rate of black hole. 



4 Asymptotic behaviors 

We consider asymptotic behaviors of our solution on the outer horizon and at infinity. 
Then, we write asymptotic amplitudes by using our solution explicitly following the 
definition given in the book of Chandrasekhar[10]. The general properties derived by 
using the asymptotic behaviors before are compared to the explicit forms derived here. 
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By doing this comparison, we can examine the structure of the solution and get a deep 
incite of our solution. The result presented here will be useful to derive various quantities 
by our solution. 



(a) Definition of asymptotic amplitudes R^ mc \ -R^ ref) and R^ ran ' ] 
The asymptotic behaviors are defined by[10] 

R s -> Rf lc) -e~ lwr * + R { I ci) -?—e iujr * (r -> oo), 

-> R(*™s)A- s e- ikr * (r->r+), (37) 

where fc is defined in Eq.(4) and r* is defined by dr*/dr = (r 2 + a 2 )/ A. Then, we find 

cur* — > z + e\nz (r — > oo), 
Ax* — > e + ln(— x) (r — > r+). 

(38) 

(b) Asymptotic amplitudes expressed by using our solution R\ n . s 
Here we express 

^(inc^ ^(ref) and ^(trans) explicitly by 

using the analytic solution. In 
below, we consider the cases of s = 0, ±1/2 (neutrino), s = ±1 (photon), ±3/2 (massless 
spin 3/2 particle) and ±2 (graviton). 

The asymptotic amplitude on the outer horizon _R( trans ) is given from Eq.(12) by 

/ fK \2s oo 

itf ans > = A s - £ <( s ) • ( 39 ) 

^ ^ / n=-oo 

To derive asymptotic amplitudes at infinity i?( mc ) and R^ ei \ we need some nontrivial 
works. We define 

R v C s -> ^. s z- 1 e-' (z+elnz) ± A!l. s z _1_2s e* (2:+eln2:) . (40) 
Then, from Eq.(17), we find 

A v n-l+s-ie i(n/2)(u+l-s) ~7ve/2 T ( U + 1 ~ S + le) ~ 

+;s " r(, + i + s -, e )i^ a « (s) ' 

a* , = 2 - i - s+ ^e-^/ 2 ^ +i+s ^e-^ 2 y (-ir ^ + 1 + s ~' e ? w <( g ) . 



(41) 
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Since the asymptotic behaviors of R\ n are expressed by the combination of Rq s and 
i?^" 1 ) the asymptotic amplitudes look quite involved. It is quite important that these 
asymptotic amplitudes are expressed in simpler forms to examine the T-S identities by 
using the asymptotic behaviors, to derive the absorption rate and the numerical compu- 
tations. We found various useful relations. Firstly, we relate A v ,. s to A+Yf 1 . By using 



Eq.(31), we find 



A -v-l 

AZ 1 . 



.^_^ sm7r(t/-g + ie) 



sin7r(z/ + s — it) 



ie mv A v _. s . 



(42) 



By combining Eqs.(25), (40) and (42), we find 



R (inc) 



A s e ieK 



l+2s 



sm7r(f + s — ie) 



(43) 



These are expressed in compact forms. In addition, we find the relation between A±. s to 
A u ± . s by using Eq.(32) 



A -;s 



T(v+ 1-s + ie) 



r(z/+ 1 +s + ie) 



J2n=-op a n( S ) 



2~ 2Sg — 



(44) 



These relations are useful to discuss the T-S identities. 



5 The Teukolsky-Starobinsky identities 

The Teukolsky-Starobinsky (T-S) identities [6], [7] for s = 1/2, 1, 3/2 and 2 can be ex- 
pressed as 

A s {D r ) 2s A s R s = C*R- S , (T-S identity (A)) (45) 
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(D) 2s R_ s = C S R S , (T — S identity (B)) (46) 



where D and are differential operators 



n 9 K 



uj ( d e + e_ 

+ l€K — I h I 



€K \dx X 1 — X ) ' 

= (47) 
\dz z — en z J 

* - | + 4- < 48 > 

where K = (r 2 + a 2 )uo — am. Note that K used here corresponds to — K defined in the 
book of Chandrasekhar[6]. In the above expression, (D^) 2s means that the differential 
operation applies 2s times and the same for (D) 2s . 
The Starobinsky constants [7] Cg are given as follows: 

| C 2 | 2 = (Q 2 2 + Aaujm - Aa 2 uj 2 ) [(Q 2 - 2) 2 + 36acum - 36a 2 cu 2 ] 

+ (2Q 2 - l)(96aV - A8aujm) + \Uu 2 (M 2 - a 2 ) , (49) 



(C f ) 2 = (g | -^) 2 (g | + i)-16aV(g f -^) + 16am^(g f -^), (50) 



(d) 2 = Q\- Aa 2 uj 2 + 4amuj , (51) 
{C,f = Q, + \ , (52) 

where 

Q s = E(s) + a 2 J 2 — 2aujm . 

Trivially C = 1. The coefficient C 2 is a complex number, but others can be real. The 
coefficients C 2 and C 1 are first obtained by Teukolsky and Press [7] and C| is given by 
Torres del Castillo[ll]. 
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Solutions of Teukolsky equation should satisfy the T-S identities. It is noted that if 
a function satisfies the T-S relation (A), it also satisfies the T-S relation (B), and vice 
versa. This is shown by identity 

A s {D*) 2s A s (D) 2s R_ s =\ C s | 2 R_ s , (53) 

for s = 0, 1/2, 1, 3/2, 2. The proof for s = 1 and 2 were given in the book of Chan- 
drasekhar[10] and the other cases can be proved similarly. This fact is important for our 
solution because it may be hard to prove analytically that our solution satisfies both of 
the T-S identities, although it is a solution of the Teukolsky equation. This is because 
our solution is expressed as series. 

Fortunately, we can show analytically that it^ n;±s 's for s — 0, 1/2, 1, 3/2, 2 satisfy the 
T-S identity (A). By using this, we can fix the relative normalization between solutions 
with spin weight s and — s analytically. 

Firstly, we consider the hypergeometric type expression for R\ n .± s in Eq.(12). As we 
see in Appendix A, we find 

T(u+l + s + te) 2 
T(v+l-s + ie) in; ~ s ' 

(54) 

Then, it is clear that by using the normalization factors in Eq.(13), the righthand side 
reduces to C*R\ irr _ s . That is, by taking the relative normalization in Eq.(13), the hyper- 
geometric type expression for R V U1 .± S satisfies the T-S identity (A). 

In addition, we can also show that the Coulomb type expression for R\ n .± s satisfies 
the T-S relation (A) as well. To prove this, we first observe that the following relation is 
satisfied: 

A s {D ^s A s R „ = J_ R , (55) 
UJ ZS 

for s = 0, 1/2, 1, 3/2, 2. The proof of this relation is given in Appendix A. Then, we 
find from Eq.(20) 

A S (D^A S R» C , S = -^R^ s . (56) 
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A s (D r ) 



2.s 



A, 



6K\ 
CO J 



2s 



r(i 



2te + ) 



r(l + s - 2ie + ) 



By using the Coulomb type expression in Eq.(25), we can immediately find 

A p ieK / \ 

A s {D ^s A s^ s = _^_(^ (a ) i% ._ - + jr _ v _ l(a)i25r -l) 

= C* s Rl ; _ s , (57) 

where we used the relation among K v {s) given in Eq.(35) and the relative normalization 
in Eq.(13) to derive the last equality. In summary, both expressions of R\ n . s satisfies the 
T-S identity (A). 

It is not possible to prove analytically that our solution R\ n . s satisfies the T-S identity 
(B), although it should satisfy by the identity in Eq.(53). It is expected that the T-S 
identity (B) will give some nontrivial identities if it applies to our solution. 

The analysis by using the asymptotic behaviors of the solution has been made by 
Teukolsky and Press [7]. By substituting the asymptotic behaviors in Eq.(37) on the 
outer horizon into the the T-S identity (B), one finds for s = 0, 1/2, 1,3/2,2 

C s \ u) r(l-s-2ie+) V ; 

Similarly, by considering the T-S identity (A) and (B) at infinity, one finds 

p(inc) -ins ( 2u} ) 2S p(inc) 

e £t ^ — -5 5 

p(ref) _^ -ins p(ref) 

(2cu) 2s s ' 



(59) 



The relations in Eq.(58) and (59) are essentially the same as those by Teukolsky and 
Press [7]. By requiring that the asymptotic amplitudes in Eqs.(39) and (43) satisfy the 
Teukolsky-Press relations, we can derive some identities between the sums of coefficients 
which we shall see in Sec. 7. 

Finally by using the relation in Eq.(55), we can show that the upgoing solution with 
the relative normalization in Eq.(27) satisfies the T-S identity (A). That is, we find 

A s (D^) 2s A s R" up . s = B s A s (D^) 2s A s R u _. s 

= <?;Kp;- S - (60) 
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Thus, the relative normalization between the solution with spin weight s and — s is 
determined analytically for the upgoing solution as given in Eq.(27). 



6 The conserved current and the energy conserva- 
tion 

In below, we take s to be 0, 1/2, 1, 3/2, 2. The Wronskian gives the conserved current[7] 
which is written by 



dr* 



dr* 



J r=r+ 



dr\ 



(61) 



where Y(s) = A s / 2 (r 2 + a 2 ) 1 / 2 it!i n;s . By substituting the asymptotic behaviors in Eq.(37), 
one finds 



(R ( : ac) )(R ( i n s c) y = (R<? ct) )(R-?y - ^ s + T +) (R ( ^ ns) )(R ( T s) y 



By using the Teukolsky-Press relations in Eqs.(58) and (59), we obtain 

(2uj) 4s 



\R<s 



inc) I 2 



\C. 



12 I s 



* ei) \ 2 + S s \Rf™ ns) \ 2 , 



where 



5, = -ie ms ^ s - 2 — 



T(s + 2ie + ) 



Explicitly, we obtain 



So = 



Si = 



Ss 

2 



s? = 



1 . 



e+(e/c) 
4a/ 



(J + 4 4)(«0 2 ' 
e + (l + 4e 2 )(e«) 3 ' 



(62) 



(63) 



(64) 



(65) 
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The relation in Eq.(63) leads to the energy conservation[6],[10] 

W J-'mc;s u - C/ ref;s u -^transjs / ap\ 

dtdn = dtdn dtdn ' 1 j 

where d 2 E- mc - s / dtdVL is the incident energy going into the black hole, d 2 E rc f. s /dtdQ is the 
energy reflected by the black hole and d 2 E tra _ ns . s / dtdVt is the energy absorbed by the black 
hole. It is clear that 5 s 's are proportional to e + for bosons and are positive definite for 
fermions. Since 2e + = e + r = e(l + -) — j^, e + can be negative for large angular 
momentum of black hole. That is, the super radiance occurs for bosons for a > e ^^ M ), 
but not fermions [12]. 

By substituting the explicit forms of asymptotic amplitudes in Eqs.(39) and (43), we 
have an identity involving the sums of coefficients which we shall see in the next section. 
In Sec. 8, we shall derive a simple expression for the absorption rate. 



7 Identities involving the sums of coefficients 

The requirement that the asymptotic amplitude on the outer horizon in Eq.(39) satisfies 
the Teukolsky-Press relation in Eq.(58) leads to the identity for s = 0, 1/2, 1, 3/2, 2 

2 rx~i 



E <(-*) = \c s 



F(v+ 1 - s + ie) 



T(u+l + s + ie) 



E <(*) • (I - 1) 



(67) 



The requirement that the asymptotic amplitudes at infinity in Eq.(43) satisfy the Teukolsky- 
Press relations in Eq.(59) is automatically satisfied. This can be seen by using the relation 
among K v (s) in Eq.(35). No more identity arises from the T-S identities. 

Next, we consider the constraint from the Wronskian. By substituting the asymptotic 
amplitudes in Eqs.(39) and (43) into Eq.(63) and by using the identity (1-1), we find 



Eoo 
n= 



.<(*) 



= \c x 



r(z/+ 1- s + ie) 



r(z/+ 1 + s + ie) 



7-2) 



(68) 
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Since the derivation is lengthy, we shall give a brief proof in Appendix A. The identity 
(1-2) can be written in a different form by using the identity (1-1) as 

(oo \ / oo \ * 

E <(«) E <(-«) • (i -2)' 
n=— oo / \n=— oo / 

(69) 

The identities (1-1) and (1-2) have simple structures and will serve to make various 
formulas into simpler ones. Since our solution is given by the sums of either hypergeo- 
metric or Coulomb wave functions, various formulas given by our solution are expressed 
in complicated forms involving the sums of coefficients. Therefore, these identities are 
powerful to compute various quantities analytically and also numerically. 

Another important feature of these identities lies in the fact that they can be used as 
a good tool to test the accuracy of the computation. In order to obtain the Teukolsky 
function, we have to solve the three term recurrence relation by the e expansion or 
the numerical calculation. In the computation, the estimation of the accuracy becomes 
an important issue. The accuracy is directly tested by looking how accurately these 
identities are satisfied because the sums appearing in these identities correspond to the 
asymptotic amplitudes themselves. As for the accuracy test of the computation, to check 
the f independence of K v {s) will be another good test. These identities and relation give 
nontrivial test of the accuracy of computation. 



oo 



-oo 



(v + l-s + ie) n 11 



8 The absorption rate and the evaporation rate of 
black hole 



In this section, we present the absorption rate formula expressed by using our solution. 
Since the calculation is lengthty, we give the brief derivation in Appendix A and only 
give the result 



^j>(trans) 



-R(inc) 



(2a; 



lis 



|C S | 2 



R 



(inc) 
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(2e/c) 



2i/+l* 



3 « f sinh27re-| 
71 cosh27re_ 



Y _|_ ; (2eK) 2 '"+ 1 e' 7r ^(-l) 2a sin7r(iy-ir) / siri7r(^-£-ie) \ 2 

7T \ sin27n/ j s 



2 ' 



(70) 



where the upper column (sinh 2ire + ) is taken for bosons (s = 0, 1, 2) and the lower column 
(cosh27re + ) is taken for fermions (s = 1/2,3/2). Here, 



F(u + 1 - ir)T(v + 1 - s + ie)r(i/ + 1 + s + ie) 



r(2z/ + i)r(2i/ + 2) 



(71) 



where 



o 

E 



(-1)' 



(v + l + s-ie) 



n=—co 
oo 



X 



E 

n=0 



(-ra)!(2i/ + 2) n (i/ + l-s + ie) 
(2i/ + l) n (u + 1 + s - ie) n „ 



n v ( \ 



<(*) 



To derive the above formula, we used sin7r(z/ + s — ie) 



1 



\2.s 



sm7r z/ 



(72) 
ie). It 



may be worthwhile to note that d v s starts from 1 in the e expansion and is the only part 
which we have make a calculation by solving the three term recurrence relation. Since the 
mathematical structure is simple, it is easy to estimate the rate either in the expansion 
of e or the direct numerical calculation. 

The evaporation rate from the black hole is given by [13] 

r. 



<N> 



{ (sinh 2-7^ 
21 



i-i 



(73) 



[ (cosh2/re + ) 1 

Here, the upper column for bosons and the lower column for fermions, k' is surface gravity, 
Q + is angular velocity on horizon and we used 

2Mu - m± 



2e + = e + r = 2Muj + 



' M 



1 " (*)' 



uj — mVL A 



(74) 
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By substituting T s in Eq.(70) into the evaporation rate formula, we find 

<N> = V^i^~ mn +)(2e/^ +1 

7T 

x- 



1 ■ (2eK) 2lJ + 1 e i7n/ (-l) 2s smir(v-iT) ( sin TrQ-s-ie) \ 2 r-j 
i+? 7r V sin27rz/ J U s 



2 • 



(75) 



It is interesting to observe that the denominator in the evaporation < N > which shows 
the difference between the boson and the fermion emissions is canceled by the sinh 2e + or 
cosh2e + factor in T s . As a result, the evaporation rate shows the Boltzman like behavior 
both for both boson and fermion emissions. 

In below, we made the e expansion for T s . In Ref.l, we obtained the renormalized 
angular momentum v and coefficients up to the order e 2 . In the approximation neglecting 
the 0(e 2 ) terms, we can immediately derive the absorption rate because v — I + 0(e 2 ). 
We find 

f T -^^ U l k _ 1 (k 2 + r 2 ) 

I ^tir^nUMA: 2 + r 2 ) 

X v 2/!(2/ + l)! ) ds ' {7b) 

where the upper cases for massless scalar, photons and gravitons, the lower case for 
neutrinos and massless spin 3/2 particles and 

< = 1 + m ^(^-(7TTF) f - < 77 > 

For massless scalar case (s — 0), it is understood that n^. =1 (A; 2 + r 2 ) = 1 for / = 0. In 
the zeroth order of e, it is easy to see that this expression reproduces the formulas first 
given by Page[14] for s = 0, 1/2, 1 and 2. 



9 Summary 

To obtain the Teukolsky functions in a desired accuracy is a quite important issue for the 
gravitational wave astrophysics. The analytic solution which we reported in our previous 
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paper is expected to serve for this purpose. Since the solution is given as series of special 
functions, the examination of the property of the solution is needed. 

In this paper, we used new Coulomb type solutions and constructed the analytic so- 
lutions valid in the entire region of r for all finite e. By using this solution, we examined 
various properties of solutions of Teukolsky equation. We found that the relative nor- 
malization between the solution with a spin weight s and — s is determined analytically 
by using the T-S identities. In addition, we found an identity (1-1) involving the sums 
of coefficients of series. By using the Wronskian, we found another identity (1-2). The 
sums of coefficients appeared in these identities are proportional to the asymptotic am- 
plitudes on the horizon and at infinity and thus can be used to test the accuracy for 
the computation. The remarkable fact is that we can compute the Teukolsky functions 
systematically in a desired accuracy. By using these identities, we demonstrated that 
the asymptotic amplitudes can be expressed in compact forms. The fact that various 
quantities can be written simpler forms by using these identities is important to see the 
general properties of these quantities analytically without any approximation and also to 
make the numerical computation. 

As an example, we derived the absorption rate and the evaporation rate of black 
hole. We found that they are expressed in compact forms from which we can extract the 
physics transparently. 

Now, the properties of solutions are well understood and it is ready to go to derive 
the Teukolsky functions by either e expansion or the numerical computation which we 
are now pursuing. 
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Appendix A: Proofs of various formulas 

In this appendix, we give brief proofs of various formulas given in the text. 

(a) The Teukolsky-Starobinsky identities 
Proof of equation (54) 

For s = 0, 1/2, 1, 3/2, 2, we rewrite R\ n . s as 



oo 



E( n . s = A s e leKX {-x)- s - le+ {\ -x)- s - le - I (1 -x) s+2ie - <(*) 

I n=— oo 

xF(n + v + 1 — ir, — n — z/ — ir; 1 — s — 2ie + ; x)} . 



Then, one finds 



AS(jD t)2 WjR . i;s 

= 4,(--J e ! »H s ^(l-x) s --|(l-x) s *- x: « 
x F(n + z/ + 1 — ir, —n — v — it; 1 — s — 2ie + ; x)} . 



Here we use the relation among hyp ergeo metric function 
d 



dx n 



(l-x) a+b ~ c F(a,b; c;x) 
= (c ~ a)w(c ~ &) " (l - x) a+b - c - n F(a, b;c + n; x) 



and then the right hand side of Eq.(A.3) becomes 
A s f^-^ 3 e ieKX {-x) s - le+ {\ -x) ie - 



{—n -v - s - ie) 2s {n + v + 1 - s - ie) 2s 



x E <( s ) 



n=— oo 



1 - s - 2ie) 2s 



xF(n + v + 1 — it, —n — v — ir; 1 + s — 2ie + ; x) . 
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By using 



(-n 



s - ie)2s{n + v + 1 - s — ie) 



2s 



("I) 



2,s 



T(u + l + s + ie) 



r(i/ + l -s + ie) 



{v + 1 + s + ie)„ 



(z/ + 1 - s + ie)„ 



and the relation between a^(s) and a v n {— s) in Eq.(32), Eq.(54) is proven . 

Proof of Eq. (55) 

The outline of the proof is as follows: 

F(u+l + s-ie) K ' + ' s 

CO 

= z s (z - 6K) s (D^) 2s e- iz z- ie -(z - €K)- ie + 

n=— oo 

xz u+s+ie (2z) n V(n + v + 1 - s + ie, 2n + 2v + 2; 2is) 



72s 



x — ^ +s+ie (2,z) n *(n + v + 1 - s + ie, 2n + 2i/ + 2; 2iz) 

oo 

= u 2s e- lz z v+ie+ (z - tK) s ~ ie+ ]T *X(s) (n + v + l- s + ie) 2s 

n=— oo 

x (-n -v - s + ie) 2s (2z) n ^(n + v + 1 + s + ie, 2n + 2v + 2; 2iz) , 



where we used the relation of confluent hypergeometric functions 



(x a+n_1 ^(a,c;x)) = (a) n (a-c+ l) n x a ~ l ^(a + n,c;x) 



(A.5) 



(A.6) 



(A.7) 



to derive the last line. By a careful calculation by using Eq.(A.6) and the relation between 
a y n (s) and a v n (— s) in Eq.(32), we can derive the relation for R v + in Eq.(55). 
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Similarly, we find 



u 2s 2~' / e- 7Te e- i7T( - u+1+s) A s (D r ) 2s A s R u _. 



= e~ lz z s - le - 



[z-eK)-** f; f<( s )- (i/ + 1 + s lt)n 



x 



j2s 



dz 2s 



2iz v+s+ie 



{v + 1 - s + ie) n 
(2z) n ^{n + v + 1 + s - ie, 2n + 2v + 2; -2iz) 



x (2^) n ^(n + i/ + 1 - s - ie, 2n + 2v + 2; -2iz) , 



where we used the relation 

- (e- x x c - a+n - l ^(a )C -x 



dx r ' 



-l) n e- x x c - a -^{a-n,c] x). 



(A.8) 



(A.9) 



By using the relation between a^(s) and a v n {— s) in Eq.(32), we can derive the relation 
for R"_. s in Eq.(55). 

(b) The derivation of the absorption rate 
We find from Eq.(24) with f = 



2 7T (2eK)- 2v+2s r(l-s + 2ie 4 

1 22 I sin7r(s + 2ie + ) T(s + 2ie+) 



T(u + 1 + s + ie) 



T(u + 1 - s + ie) 



1 



(A.10) 



where .D^ is defined in Eq.(71). Then, by changing v to — v— 1 and using x (s) = a^(s), 
we find 



|K_„-i(s)| 2 (2eK)^ 



\Ku(s) 



IT' 



sin TT(V + IT) 



I sin n(u + s + ie) 
\ sin 2tcu 



(D 



s) 2 



(A.ll) 
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We also find by using Eqs.(35) and (36) 
K v {s)K- v -i{sY _ K u (s)K_ u _ 1 ( 



\K»{s) 



{2en) 



2u+l 



Sin 7T V 



' K-u-i(s) 



it) 



7Y 



s + ie) 



sin 2ixv 



(A.12) 



By using Eqs.(A.ll) and (A.12), we obtain 



K„(s)-ie— ™*<y-° + *) 
sin ix\y + s — ie) 

1 ir {2eK)- 2v+2s T(l-s + 2ie 



K-^s) 
I> 



D»2 2s simr(s + 2ie + ) T(s + 2ie H 



s + ie) 



T(v + 1 - s + ie) 



x 



i + i(-l) 2s e^(2e/t) 2 ^ +1 sin7r(z/-ir) /sin7r(i/ - s + ie) y ^ v 



71 



sin27rz/ 



(A.13) 



By using this formula, we can derive the formula for the absorption rate in Eq.(70). 

(c) Identity (1-2) form the Wronskian 

Here we show that the identity (1-2) is obtained from the Wronskian given in Eq.(63). 
Let us first parameterize |i?^ mc ^| 2 as 

\A s \ 2 2 2s - 2 e-*' 



\R { } nc) \ 2 



-{{K^+ie^K^s) 



+ 



K v (s) - ie 



,sin7r(z/ — s — ie 



riJU-i(s) 



\K v {s)+ij™K- v _ 1 {s)\ 2 } 



sin7r(z/ + s + ie) 

T(u + 1 - s + ie 



T(u + 1 + s + ie) 



n=— oo 



(A.14) 



Then, by using the relation 

sin tt(v — s — ie 



\K v (s) - ie 



sin7r(z/ + s + ie) 



r^-K^is)] 2 - \K v {s)+ie^K_ v _ 1 {s) 



-2ie lns (eK) 2s+1 e m 



r(l - s + 2ie+) 
T(s + 2ie + ) 



r(z/+ 1 + s + ie) 
T(u+ 1- s + ie) 



(A.15) 
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we can prove that the contribution from the second part in the parenthesis in Eq.(A.14), 
i.e., the part in Eq.(A.15) cancels the transition (absorption) part <5 s |i?^ trans ) | 2 . Thus, 
the contribution from the first part in Eq.(A.14) should be equal the reflection part. By 
equating these two, we can derive the identity (1-2). 
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